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We present a new method of constructing the matrix product (MP) states and apply it to the 
ground state of the Heisenberg spin chain with ahernating spins 1 and i and antiferromagnetic 
exchange interactions between nearest neighbors (the simplest example of a quantum ferrimagnet). 
The elementary matrix state is constructed in a way which ensures given transformational prop- 
erties under rotations, which allows one to fix the total spin and its z-projection of the entire MP 
wavefunction. We compare the variational MP results with the numerical results obtained through 
a Quantum Monte Carlo method; the agreement is found to be within 0.4% for the ground state 
energy and 5% for the correlation functions. 
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Recently, there has been a considerable progress in 
the study of one-dimensional (Id) quantum spin chains 
by meanS|-p£. the so-called matrix product (MP) states 
techniqueJjij For the first time, they were introduced^ as. 
a convenient representation of the valence bond statesu 
(VBS) for integer-spin chains. In terms of the original 
spin states, VBS wavefunctions are highly non-local, but 
they factorize in terms of matrix states which consid- 
erably simplifies all calculations. Shortly after that, MP 
states have proved to be a useful tool for constructing new 
classes-pf|MDin Hamiltonians with exactly known ground 
states,0~El'l3 for variajtional study of the ground state prop- 
erties of spin chaina] and spin ladders ,u and also for the 
study of elerp.ea.tary excitations which are solitons in the 
string order 

In this paper, we present a new method of construct- 
ing the matrix product (MP) states with fixed quan- 
tum numbers, namely the total spin and its z-projection, 
and demonstrate its use on the example of the Heisen- 
berg spin chain with antiferromagnetic nearest neighbor 
interaction and alternating spins 1 and i This sys- 
tem, according to the Lieb-Mattis theorcmjiil has a de- 
generate multiplet of ground states with the total spin 
S'tot — L/2; here L is the number of unit cells, each cell 
consisting of spins 1 and ^ . In each of the ground states 
the rotational symmetry is spontaneously broken, and 
the long-range magnetic order exists jl3 thus, one may 
consider this model as the simplest example of a quan- 
tum ferrimagnet. It should be mentioned that there ex- 
ist real magnetic materials belonging to the family of 
Cu(II)Ni(II) complexes, which are well described by the 
above model, see Ref. and references therein. For this 
model, we construct a MP wavefunction which has only 
two variational parameters, and show that our variational 
results for the ground state properties are in a very good 
agreement with the numerical data obtained by means of 



the Quantum Monte Carlo (QMC) technique. 

Matrix states with given quantum numbers. First, let 
us recall some basic facts about the MP states. A matrix 
product state is defined as 



l^^> =Tr(5ig2---5L), 



(1) 



where the elementary matrix states gi are matrices com- 
posed from the spin states of the i-th magnetic elemen- 
tary cell. The dimension of the elementary matrix de- 
pends on the problem: for example, translationally in- 
variant VBS state for spin-5' chain can be prepresented 
with the help of (S' -I- 1) x (S' -I- 1) matrices.El In the fol- 
lowing, we will assume that gi is a square matrix. It is 
worthwhile to note that since the definition of (|l|) involves 
the trace sign, the MP state does not depend on the rep- 
resentation of the elementary matrices, i.e., any global 
unitary transformation gi i— > UgiU^ keeps \n) invariant. 

Let us demand that the elementary matrix g has cer- 
tain transformational properties. Namely, we would like 
to construct a multiplet of matrices g^"^ transforming 
according to the representation of the rotation group 
under rotations Tf? : 



Tng- 



jm 



(2) 



where the sign ~ means unitary equivalence. Further, as- 
sume that the desired dimension of the elementary ma- 
trix is {N+l)x{N+l) and that we know the complete 
set of wavefunctions of the elementary cell {IV'a^)}, A 
being the total spin and fi its z-projection. One can 
choose a basis in the matrix space consisting of the ma- 
trices X'"i = {T'"iy, k = 0...N, q = -k...k. Here 
T'^'^ are the (7V-|-l)x(iV+U-| matrix representations of 
irreducible tensor operatoraiI'LJ transforming according 
to the I?*^ representation of the rotation group; defining 
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j^kq _ j-j-feg^t jygt ensures that X'"^ transform according 
to the conjugate representation V'', Hke usual wavefunc- 
tions. In the simplest case A^=l (2x2 matrix space) the 
T^'^ matrices are just the identity matrix 1 and the Pauli 



matrices a 
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explicit expressions for T*-'^ in 3 x 3 space are given below. 
Then, our basis transforms as 



kq' 



(4) 



which on the other hand amounts to a unitary transfor- 
mation TrX'"! = U{R) X'"iW[R) and does not have any 
effect on the MP wavefunction. 

It is easy to show that the following construction sat- 
isfies our requirement (||): 



(5) 
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Here {jm\kq, Xfi) are the Clebsch-Gordan coefficients, 
and c^^ are arbitrary constants. In fact, we used the 
usual quantum mechanical addition rules for angular mo- 
menta, demanding that g has certain quantum numbers 
j, m characterizing the "hyperspin" defined in an ex- 
tended space which is a direct product of (7V+1) x (iV+1) 
matrix space and the Hilbert space of the cell. 

Then, we can construct a translational invariant ma- 
trix state with given hyperspin J = Lj and its z- 
projection M = J: 



(6) 



The matrix state G'^'^ belongs to the multiplct transform- 
ing according to V'^ and thus has the same structure (||) 
as g^'' ^ but now have the meaning of many-body 

wavefunctions of an open spin chain (their total spin 
A = |J — A^|,...,J-|-A^ can differ from J because of the 
edge effects). Taking the trace \^) = Tr(G''^'^), we pick 
up only the component, i.e. the state \tpjj) describ- 
ing a chain with periodic boundary conditions. Other 
states of the multiplet with M < J can be obtained from 
Itpjj) by applying 5'^^. 

It should be mentioned that the construction (^ al- 
lows several different values of j even for the same size 
of the magnetic elementary cell; the number of possibil- 
ities grows fast with increasing the cell size, when more 
and more complicated magnetic structures with longer 
period are allowed. One can also construct various trans- 
lationally non-invariant 'excited' MP states, but this goes 
beyond the scope of the present paper. |^ 

One can see that, known MP states for VBS modelsEI 
and spin-i laddersQ indeed follow the general structure 
described above. For example, the elementary matrix for 
the spin-1 AKLT state can be written in the form 



<7AKLT = l/V3{(a+)t|-> + (a-)t|+) - (a«)t|0)} , 

which is exactly g'^" in our notation, and four components 
of the matrix defined by (^) describe a singlet and 
the Kennedy triplet states of the open spin-1 chain. 

In the spin-2 case one possibility is to construct T^™ 
as the matrix representations of S=l spin operators, and 
T^™ can be expressedllS as bilinear forms of T^"\ This 
yields the following expressions for T'"^ in 3x3 space: 




(7) 



The matrix T^^ is proportional to the identity matrix 1, 
and the relationship T^'-™ = (-l)™(rJ™)t holds. If the 
elementary cell consists of a single 5=2 spin, then the 
elementary matrix g^'^ constructed according to (^ is 
5S=2 = E™ff^")^lV'2m), which coincides with the 5=2 
VBS matrixp nine components of describe degener- 
ate ground states of the open spin-2 VBS chain. 

In case of a spin-^ ladder the elementary cell consists 
of two spins, and according to (||) one gets the ansatz 
used in Ref. 0: 



^/ladder 



„00 



1 . |,S) 



C^SAKLT • 



Here \s) denotes the singlet state, and the states \^) in 
.9AKLT, /i = 0, ±1 now should be interpreted as the triplet 
states of two S=^ spins. 

Our ansatz (||) allows one to construct MP states with 
spontaneously broken rotational symmetry (nonzero to- 
tal spin J), which is practically important for the case 
of ferrimagncts. In the next Section we apply our for- 
malism to the simplest model of a quantum ferrimagnet 
and argue that it really provides a good description of its 
ground state properties. 

Application to the Heisenberg model with alternating 
spins 1 and i. We consider the model described by the 
Hamiltonian 



(8) 



where Sn and t„ are respectively spin-1 and spin-^ oper- 
ators at the n-th elementary magnetic cell (with eigen- 
states denoted as (-|-,0, — ) and (t,i), respectively). The 
elementary cell consists of spins 1 and ^, the exchange 
interaction exists only between nearest neighbors, and 
all bonds are antiferromagnetic and of the same strength 
which is set to unity (see Fig. Q). 
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According to the Lieb-Mattis theorem|l3 the ground 
state of this model has total spin ^tot = L/2, where L is 
the number of elementary cells. The rotational symmetry 
is spontaneously broken in the ground state, which leads 
to a ferrimagnetic long-range order. Thus, we have to 
construct the elementary matrix g with j = 1/2, ni = 
1/2. The complete set of cell wavefunctions il^x^ contains 
one doublet (A = 1/2) and one quartet (A = 3/2): 



n n+1 

FIG. ^. The Heisenberg spin chain with alternating spins 1 
(large circles) and i (small circles); boxes show the unit cells. 



IV'l/2.1/2)-V^I+i>-l/\^|0T), 

IV'i/2,-1/2) = T) + l/\/3 |0 i) , 

IV'3/2,3/2) = 1+ T) , 1^3/2,-3/2) = I- i), (9) 
IV'3/2,1/2) = V273|0T) + l/V3|+i), 

IV'3/2,-1/2) - V273I0 i) + i/Vi\- T) , 

Let us choose the simplest case of 2x2 matrix space, 
1, = a", Xi^±i - -(jT (cf. 
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then the basis is X 

Eq. (^), and, according to (||), the most general form of 
the elementary matrix with j ~ m = 1/2 is 

- 1 1^1/2,1/2) , (10) 
= _i/^^o 1^^/2,1/2) - VV^a-' IV'1/2,-1/2) , 
= -i/Vs^o 1^3/2,1/2) - l/V6a-' lV'3/2,-1/2) 
-l/x/2a+i 1^3/2,3/2) ■ 

Thus, g contains two independent variational parameters 
u = C1/C3 and V = C2/C3 which are assumed to be the 
same throughout the chain. By the construction, the 
variational MP state |il) = Tr (51(72 • ■ ■ 5l) has total spin 
'S'tot = L/2, i.e., the same as the true ground state. 

The quantum averages can be calculated using the 
transfer matrix technique; for example, for any opera- 
tor F acting on the states of only one unit cell one has 
{n\F\n) = Tr(G^-iGF)/Tr(G^), where G = g* ^ g 

is the transfer matrix, Gp — g* ^ {^9)1 ^-nd the sign 
(g) means an outer matrix product; further details can 
be found in Refs. The resulting expression for the 
ground state energy per unit cell is rather lengthy but 
still manageable. 
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£;var = {1/18Q){A + B/Z), 

Z=[Uv\u^ + l) + 3]'/^, Q 

A = -180m'* - SGOu^v'^ - 20v'^ 

~180V2u^v + SGVGuv^ - 40\/2u^ 

-108m^ + SeVSuv ~ 96i;2 -I- 27V6u - 29V2v 



B = - 1440m V -480w2u'* 



-1081;"- 210 V2u\ 



180V2u'^v 



h 108\/6u"i; 
288V3uw^ 



(11) 



4, 

3„.2 



234V6WW 



'144u2 -I- lOSVSuv - lOSv'^ 



72V2v^ 
81a/6u - 27\/2u -I- 18 . 



Its numerical minimization gives -Eva'" = ^mp = —1.449 
(at u ~ —1.303, V — 1.079). Using this solution, one can 
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FIG. 1^. Longitudinal spin correlation functions K^^t (I), where 
s and s' can be either S=l or S=^, and / is the distance 
between spins in unit cell lengths, so that for integer (half- 
integer) / the correlation is measured between spins of same 
(different) type, respectively; (a) correlations between 5=1 
and any other spin; (b) correlations between 5=| and any 
other spin. The dashed line shows the asymptotic value pre- 
dicted by the spin wave theoryJlS The QMC data were ob- 
tained at L = 32. The numerical uncertainty is less than 
0.002, which is much smaller than the symbol size. 

also calculate average spin values (sf ) and spin correla- 
tion functions (sfs-^„), where s,s' € {5, r}. Correlation 
functions behave as a -I- 6e^"/^, the correlation length ^ 
being extremely short, ^ = 0.365 in the unit cell lengths. 
In the regular spin wave (SW) theory the decay is ex- 
ponential only asymptotically, and ^ « 0.7|i3 the spin 
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reduction is also overestimated in the SW calculation. 

The results were compared with the numerical data ob- 
tained with the quantum Monte Carlo (QMC) method 
based on the Trotter-Suzuki decomposition of checker- 
board type. The raw data for a set of the Trotter num- 
bers n were extrapolated into the n oo limit with a 
parabolic fitting formula. We have calculated the chains 
of length L = 4, 8, 16, and 32 under the periodic bound- 
ary condition in the subspace with fixed total magne- 
tization = L/2. Though the L = 32 data are al- 
ready close to the corresponding thermodynamic limit 
values, we have carried out the L — > oo extrapolation. 
We have checked that almost the same results are ob- 
tained at two temperatures, T — 0.04 and 0.02, and 
thus we assume that T — 0.02 is low enough to de- 
scribe the ground state properties successfully. Typically, 
about 10^ MC steps were needed to reach the equilib- 
rium. The ground state energy per unit cell was esti- 
mated to be i?QMC = —1.455 ± 0.001 in the thermody- 
namic limit, which agrees very well with the MP varia- 
tional result. The agreement is also good for the average 
spin values, the MP calculation gives {S^} = 0.779 and 
(rf) = —0.279, to be compared with the QMC results 
0.793 ± 0.002 and -0.293 ± 0.002, respectively. In Fig. | 
we show the results for various spin-spin correlation func- 
tions obtained in the MP approach, together with the 
QMC data. The discrepancy never exceeds 5%,|-S¥hich is 
much better than the corresponding SW result .113 

Summary. To conclude, we propose a new version of 
the matrix product (MP) states approach to the descrip- 
tion of quantum spin chains, which allows one to con- 
struct MP states with certain quantum numbers J (the 
total spin) and M (its ^-projection). Known MP rep- 
resentations fof VBS states for integer-spin antiferromag- 
netic chainSjtl and recently proposed MP ansatz for S=^ 
laddeiH correspond to the particular case J = 0, M = 
of our ansatz. The method can be useful for the varia- 
tional description of quantum ferrimagnetic chains whose 
ground state has nonzero J, and of higher-S* spin lad- 
ders. We apply our approach to the isotropic ferrimag- 
netic Heisenberg chain with alternating spins 1 and ^ ; on 
the other hand, we study this system numerically using 
the quantum Monte Carlo (QMC) technique. For both 
the ground state energy and the correlation functions we 
obtain very good agreement between the variational re- 
sults and the QMC data. We think that the fundamental 
reason for the success of the MP approach to the ferri- 
magnetic chain is that in this model the correlations of 
fluctuations above the long-range order are very short- 
ranged. MP ansatzes generally yield short range correla- 
tions which is just appropriate here. 
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FIG. 1. The Heisenberg spin chain with alternating spins 
1 (large circles) and | (small circles); boxes show the unit 
cells. 

FIG. 2. Longitudinal spin correlation functions K^^,{1), 
where s and s' can be either S=l or S=^, and / is the dis- 
tance between spins in unit cell lengths, so that for integer 
(half-integer) I the correlation is measured between spins of 
same (different) type, respectively; (a) correlations between 
S=l and any other spin; (b) correlations between and 
any other spin. The dashed line shows the asymptotic value 
predicted by the spin wave theoryJlj The QMC data were 
obtained at L = 32. The numerical uncertainty is less than 
0.002, which is much smaller than the symbol size. 
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